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We extend the confined covariant constituent quark model that was previously developed by us for 
mesons to the baryon sector. In our numerical calculation we use the same values for the constituent 
quark masses and the infrared cutoff as have been previously used in the meson sector. In a first 
application we describe the static properties of the proton and neutron, and the A-hyperon (magnetic 
moments and charge radii) and the behavior of the nucleon form factors at low momentum transfers. 
We discuss in some detail the conservation of gauge invariance of the electromagnetic transition 
matrix elements in the presence of a nonlocal coupling of the baryons to the three constituent quark 
fields. 
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I. INTRODUCTION 

We use the confined covariant constituent quark model (for short: covariant quark model) as dynamical input to 
calculate the electromagnetic transition matrix elements between light (u, d, s) baryons. In the covariant quark model 
the current-induced transitions between baryons are calculated from two-loop Feynman diagrams with free quark 
propagators in which the high energy behavior of the loop integrations is tempered by Gaussian vertex functions [iH^ . 
An attractive new feature has recently been added to the covariant quark model inasmuch as quark confinement has 
now been incorporated in an effective way, i.e. there are no quark thresholds and thus no free quarks in the relevant 
Feynman diagrams ^ . We emphasize that the covariant quark model described here is a truly frame-independent 
field theoretical quark model in contrast to other constituent quark models which are basically quantum mechanical 
with built-in relativistic elements. 

In the covariant quark model we use the same values for the constituent quark masses and the infrared cutoff for 
all hadrons (mesons and baryons) independent of the hadron masses. We believe that the formulation of the confined 
covariant quark model constitutes a major advance both from the conceptual and the practical point of view. While 
an unconfined quark model is valid only for hadrons with masses low enough to lie below the sum of the constituent 
quark masses the confined covariant quark model can be applied to all hadrons regardless of their masses. The 
viability of the improved covariant quark model was demonstrated in a number of applications to mesonic transitions 
in 0]. The form factors of the B{Bs) P{V) transitions were evaluated in in a parameter- free way, in the full 
kinematical region of momentum transfer. As an application, the widths of some nonleptonic Bg decays were also 
calculated. This approach was successfully applied to a study of the tetraquark state X(3872) and its strong and 
radiative decays (see, Refs. 0, @|)- 

In the present paper we formulate the covariant quark model with infrared confinement for the baryon sector. By 
keeping the same values for the constituent quark masses and the infrared cutoff as had been used in the meson 
sector we are able to reduce the number of free model parameters in the baryon sector to essentially the set of baryon 
size parameters. As a first application we describe the static properties of the nucleon and the A-baryon (magnetic 
moments and charge radii) and the behavior of the nucleon form factors at low momentum transfer. In a forthcoming 
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publication [9| we are planning to study the rare decays of the Af,-baryon into the or the neutron. The present 
paper provides the necessary input for such a calculation in as much we determine the properties of the light baryons 
from their electromagnetic interactions. 

The paper is structured as follows. In Sec. II we review the basic notions of our dynamical approach — the covariant 
quark model for baryons. We present the interaction Lagrangian describing the nonlocal coupling of a proton to its 
constituents, discuss the choice of interpolating currents and the vertex function, recall the compositeness condition 
for bound-state hadrons and show how the confinement ansatz is implemented in the baryon sector. In Sec. Ill we 
include the electromagnetic interactions of quarks and charged baryons in a manifestly gauge-invariant way, derive 
the Lagrangian describing the nonlocal interaction of the baryon, quark and electromagnetic fields. In Sec. IV we 
present the loop integration techniques that allow one to calculate the nucleon mass function and its derivative and 
the electromagnetic form factors of the nucleons. By analytically verifying the pertinent Ward and Ward-Takahashi 
identities we discuss in some detail how gauge invariance is maintained in the electromagnetic transitions. Sec. IV 
also contains our numerical results for the magnetic moments and form factors of the proton and neutron. We find 
that a particular superposition of vector and tensor interpolating currents gives satisfactory results for the nucleon 
static properties and form factors at low energies. In Sec. V we extend our approach to describe the static properties 
of the As hyperon. We summarize our findings in Sec. VI. 



II. THE COVARIANT QUARK MODEL FOR BARYONS 



The basic ingredients of the covariant quark model for baryonic three quark states prior to the implementation of 
confinement can be found in [ll-0|- This includes a description of the structure of the Gaussian vertex factor, the 
choice of interpolating baryon currents as well as the compositeness condition for baryons. 

The new features introduced to the meson sector in 0, 0] and applied to the baryon sector in this paper are 
both technical and conceptual. Instead of using Feynman parameters for the evaluation of the two-loop baryonic 
quark model Feynman diagram we now use Schwinger parameters. The technical advantage is that this leads to a 
simplification of the tensor loop integrations inasmuch as the loop momenta occurring in the quark propagators can 
be written as derivative operators. Furthermore, the use of Schwinger parameters allows one to incorporate quark 
confinement in an effective way. Details of these two new features of the covariant quark model have been described 
in 1,11. 

Let us enumerate the number of model parameters that are needed in our approach for the description of baryons. 
As stated in the introduction the values of the constituent quark masses and the universal confinement parameter 
are taken over from the meson sector. The coupling strength of a baryon to its constituent quarks is fixed by the 
compositeness condition. This leaves one with one size parameter for each baryon. For the present paper we need 
the size parameters of the proton, neutron and Ag. Naturally we use the same size parameter for the proton and the 
neutron. The size parameters are determined by a fit to the static e.m. properties and form factors of the nucleons 
and the As . We have added one more parameter to the set of basic parameters which describes the mixing between 
vector and tensor interpolating currents. 

A. Lagrangian and three-quark currents 

Let us begin our discussion with the proton. The coupling of a proton to its constituent quarks is described by the 
Lagrangian 

^Fnti^) = 9NP{x) ■ Jp{x) +gN Jp{x) -pix) , (1) 



where we make use of the same interpolating three-quark current Jp{Jp) as in Ref. 11 1 



Jp{x) = ldxildX2ldX3FN{x;Xi,X2,X3)J^f{xi,X2,X3). 



'3q 

4l'{xi,x2,x3) = rVrf"H^i)-[e°'°^"'w"n2^2)cr^u'^Ha;3)], 



(2) 



Jp{x) = ldxijdx2jdx3FN{x;Xi,X2,X3)J^^{xi,X2,X3). 



J^fixi,X2,X3) = [e'^^'^^'^^^2''Hx3)r^CrHx2)] •d'^Ha^OlT'' 
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The matrix C = 7^7^ is the usual charge conjugation matrix and the {i = 1, 2, 3) are color indices. There are two 
possible kinds of nonderivative three-quark currents: (g) = 7" <8) 7q (vector current) and ig) = i a'^^ (g) 
aa/3 (tensor current) with a°'^ = ^{'f"^^ — 7^7")- The interpolating current of the neutron and the corresponding 
Lagrangian are obtained from the proton case via p — > n and u ^ d. As will become apparent later on, one has to 
consider a general linear superposition of the vector and tensor currents according to 

Jn = xjj^ + {1 - x)J^ , N^p,n (3) 

where the mixing parameter x extends from zero to one (0 < a: < 1). When taking the nonrelativistic limit of the 
vector and tensor currents one finds that the two currents become degenerate. The limiting currents for the proton 
and the neutron read 

Jj ^ 4 = e^^^^^- ar/ iru' T2aru') , Jj ^ JX = e^^^^^^ a 0^0 , (4) 

where are the upper components of the respective Dirac quark spinor fields and where Ui are Pauli spin matrices. 

Most of the properties of the nucleons are only weakly dependent on the choice of interpolating currents. However, 
in order to get the correct value for the charge radius of the neutron, one needs to use the superposition of currents 
Eq. even though the currents J"^ and become degenerate in the nonrelativistic limit. In view of the fact 
that a nonrelativistic description of the neutron gives zero values for the charge radius of the neutron [l^ this is an 
indication that relativistic corrections play a crucial role for the desciption of the neutron charge radius. 

The vertex function i^^r characterizes the finite size of the nucleon. We assume that the vertex function is real 
and the same for the proton and the neutron. To satisfy translational invariance the function i^jv has to satisfy the 
identity 

Fjv(x + a; + a, a;2 + a, + a) — Fn{x;xi,X2,X3) (5) 
for any given 4- vector a . We use the following representation for the vertex function 

Fn{x;Xi,X2,X3) = S^'^^X-'^WiXi) ^Ni'^iXi- Xj)'^] , (6) 

1=1 i<j 

where is the correlation function of the three constituent quarks with the coordinates xi, X2, X3 and masses mi, 
TO2, respectively. The variable Wi is defined by Wi = mi/{mi + m2 + m^) such that X^^^^ = 1. Note that 
F/v(a;; cci, X2, cca) is symmetric in the coordinates Xi, i.e. symmetric under Xi o Xj. 

We shall make use of the Jacobi coordinates pi^2 and the CM coordinate x which are defined by 



Xl = 2; + 75 ^3 Pi - ^ i'2W2 + W3) P2 , 



X2 = ^ + 71 ^3 Pi + ^ (2u;i + W3) P2 , 

2^3 = S - ^ (lOl + W2) Pi + 7f {Wl - W2) P2 . (7) 

The CM coordinate is given by a; = X]i=i '^iXi- In terms of the Jacobi coordinates one obtains 

- x,)^ = pj + pI . (8) 

i<j 

Note that the choice of Jacobi coordinates is not unique. With the above choice Eq. ([7]) one readily arrives at the 
following representation for the correlation function in Eq. ([6]) 

$^(5:(x.-x,)^) = |^|^e-^^(— )-''^(— )$A.(-Pf-P|), (9) 
^Ni~P?-Pi) - I jd^p,jd^p2e^''^P^+^P^P^<^^{pl + pl), 

Even if the above choice of Jacobi coordinates was used to derive © the representation © in its general form can 
be seen to be valid for any choice of Jacobi coordinates. The particular choice ([7]) is a preferred choice since it leads 
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to the specific form of the argument —Pi — = +P2 + PiP2)- Since this expression is invariant under the 

transformations: Pi ^ P21 P2 ^ ~P2 ~Pi and pi — >■ —pi — P2, the r.h.s. in Eq. Q is invariant under permutations of 
all Xi as it should be. 

In the next step we have to specify the function $jv(— P;^ — P|) = <i>Ar(— P^), which characterizes the finite size of 
the baryons. We will choose a simple Gaussian form for the function $Ar 

$jv(-P')-exp(PVA2,), (10) 

where An is a size parameter parametrizing the distribution of quarks inside a nucleon. Note that we have used 
another definition of the A^v in our previous papers: Ajv A°j^'^ / (3^/2) . 

Since P ^ turns into — Pj in Euclidean space the form ([T(I| has the appropriate fall-off behavior in the Euclidean 
region. We emphasize that any choice for $Ar is appropriate as long as it falls off sufficiently fast in the ultraviolet 
region of Euclidean space to render the corresponding Feynman diagrams ultraviolet finite. The choice of a Gaussian 
form for $^ has obvious calculational advantages. 

The coupling constants are determined by the compositeness condition suggested by Weinberg ,12,] and Salanr |l3l| 
(for a review, see [l3|) and extensively used by us in previous papers on the covariant quark model (for details, see jl5|). 
The compositeness condition postulates that the renormalization constant of the bound-state wave function is set equal 
to zero. In the case of a baryon this implies that 

Zjv = l-I]^(TOAr) =0, (11) 

where S'^ is the on-shell derivative of the nucleon mass function E^v, i.e. S'^ = dTjpf/d^, at p^ = and where 
rriAT is the nucleon mass. The compositeness condition is the central equation of our covariant quark model. It can 
be viewed as the field theoretic equivalent of the wave function normalization condition for quantum mechanical wave 
functions. The physical meaning, the implications and corollaries of the compositeness condition have been discussed 
in some detail in our previous papers (see e.g. [51). 



B. Infrared confinement 

In we have shown how the confinement of quarks can be effectively incorporated in the covariant quark model. In 
a first step, we introduced an additional scale integration in the space of Schwinger's a-parameters with an integration 
range from zero to infinity. In a second step the scale integration was cut off at the upper limit which corresponds to 
the introduction of an infrared (IR) cutoff. In this manner all possible thresholds present in the initial quark diagram 
were removed. The cutoff parameter was taken to be the same for all physical processes. Other model parameters 
such as the constituent quark masses and size parameters were determined from a fit to experimental data. 

Let us describe the basic features of how IR confinement is implemented in our model. All physical matrix elements 
are described by Feynman diagrams written in terms of a convolution of free quark propagators and the vertex 
functions. Let n and m be the number of the propagators and vertices, respectively. For the current-induced baryon 
transitions or the derivative of the mass function one has four quark propagators and two vertex functions, i.e. one 
has n = 4 and m = 2. In Minkowski space the two- loop diagram will be represented as 

„ m n 

nbi,...,p™) = [d'k]^ n *n+n(-^'+„) n 5^3(^.3+^.3), 

•' 11 = 1 .3=1 

where the vectors ki are linear combinations of the loop momenta ki. The Vi are linear combinations of the external 
momenta pi . The strings of Dirac matrices appearing in the calculation need not concern us here since they do not 

m 

depend on the momenta. The external momenta pi are all chosen to be ingoing such that one has ^ Pi = 0- 

i=l 

Using the Schwinger representation the local quark propagator is written as 

00 

S{k) = = (m+ 1^) f dp e~P . (13) 

m— f J 


As mentioned before one takes the Gaussian form or the vertex functions, i.e. 

$,+„(-X2) =exp[A+„A'2] i = i,. (14) 
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where, as in pO)) . the parameters /3i+„ are related to the respective size parameters of the baryons via /3i+„ = 1/Af . 
The integrand in Eq. (|12p has a Gaussian form with the exponential factor {kak+2kr+R) where, in case of the baryonic 
two- loop calculation, a is a 2 x 2 matrix depending on the parameters f3i, r is a. 2-component vector composed from 
the external momenta, is a 2-component vector of the loop momenta of the form k — (fci, fc2) and i? is a quadratic 
form of the external momenta. Tensor loop integrals are calculated with the help of the differential representation 

^Mg2fcr^ 1 a ^2fc,_ (15) 

After doing the loop integration the differential operators d/drj will give cause to outer momenta tensors which, in 
the present case, are p and p'. We have written a FORM [Ifl] program that achieves the necessary commutations of 
the differential operators in a very efficient way. 

After doing the loop integrations one obtains (n denotes the number of propagators) 

oo 

n = y d"/3i^(/3i,...,/3„), (16) 



where F stands for the whole structure of a given diagram. The set of Schwinger parameters /3i can be turned into a 
simplex by introducing an additional i-integration via the identity 

°? 

1= / dt^{t-Y,P^) (17) 
^=1 

leading to 

n = Jdtt"-^Jd"a6(l-J2a,^F{tai,...,tan). (18) 


There are altogether n numerical integrations: (n — 1) a-parameter integrations and the integration over the scale 
parameter t. The very large t-region corresponds to the region where the singularities of the diagram with its local 
quark propagators start appearing. However, as described in [1], if one introduces an IR cutoff on the upper limit 
of the i-integration, all singularities vanish because the integral is now analytic for any value of the set of kinematic 
variables. We cut off the upper integration at 1/A^ and obtain 

l/A^ 1 

n"= / dtt"-W(i"a^(l-^a,) i^(toi,...,to„). (19) 



By introducing the IR cutoff one has removed all potential thresholds in the quark loop diagram, i.e. the quarks are 
never on-shell and are thus effectively confined. We mention that an explicit demonstration of the absence of a two- 
quark threshold in the case of a scalar one-loop two-point function has been given in Ref. [5]. We take the infrared 
cutoff parameter A to be the same in all physical processes. The numerical evaluations have been done by a numerical 
program written in the FORTRAN code. 



III. ELECTROMAGNETIC INTERACTIONS 



We use the standard free fermion Lagrangian for the baryon and quark fields: 

£feee(a;) = B{x){i ^ - mB)B{x) + ^ q{x){i ^ - my)q{x) , (20) 

where is the constituent quark mass. The interaction with the electromagnetic field has to be introduced both at 
the baryon and the quark level. In a first step we gauge the free Lagrangians Eq. (|20p of the quark and baryon fields 
in the standard manner by using minimal substitution: 



Qt'B (a^ - ieBA^')B , ^^'q^ ^ (9^ - ie^^A^)?, , 



(21) 
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where es is the electric charge of the baryon B and Cq- is the electric charge of the quark with flavor g^. The interaction 
of the baryon and quark fields with the e.m. field is thus specified by minimal substitution. The interaction Lagrangian 
reads 

= eBBix)4Bix)+J2eqq{x)4q{x) . (22) 

As will become apparent further on, the electromagnetic field does not directly couple to the baryon fields as a result 
of the compositeness condition. 

Next one gauges the nonlocal Lagrangian E g dp . The gauging proceeds in a way suggested in Refs. [3 and 
used before by us (see, for instance, Refs. pH Il9j|). In order to guarantee local invariance of the strong interaction 
Lagrangian one multiplies each quark field q{xi) in with a gauge field exponential. One then has 

q,{xi} ^ e-''^^'^---^PU,{x^) , (23) 

where 

Xi 

I{x,,x,P) = Jdz^Af^iz). (24) 

X 

The path P connects the end-points of the path integral. 

It is readily seen that the full Lagrangian is invariant under the transformations 

^ e-../(-)9,(a;) , q,{x) ^ Ux)e-'-nf(-) , 

B{x) e*^«^'(^) S(a;) , B{x) B{x) e-"'^^^''\ 
A^'ix) ^ A^'{x) + d^'f{x), 

3 

where cb = ^qi- 

i=l 

One then expands the gauge exponential up to the requisite power of tqA^ needed in the perturbative series. This 
will give rise to a second term in the nonlocal electromagnetic interaction Lagrangian Superficially it 

appears that the results will depend on the path P taken to connect the end-points in the path integral in Eq. (|24p . 
However, one needs to know only the derivatives of the path integral expressions when calculating the perturbative 
series. Therefore, we use the formalism suggested in [l3,[l8|] which is based on the path-independent definition of the 
derivative of I{x, y, P): 

\im dx^^Iix,y,P) ^ \im [Iix + dx,y,P')~Iix,y,P)] (25) 
where the path P' is obtained from P by shifting the end-point x by dx. The definition ((25)) leads to the key rule 



^/(a;,y,P) = A^(a;) (26) 



which in turn states that the derivative of the path integral I{x, y, P) does not depend on the path P originally used 
in the definition. 

As a result of this rule the Lagrangian describing the nonlocal interaction of the baryon, quark and electromagnetic 
fields to the first order in the electromagnetic charge reads 

^c„.-„o„ioc(^) = gsBix) ■fdyA^iy)J^_^^{x,y)+gB [ dyA^{y)J%_^^{x,y) ■ B{x) , (27) 
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where the nonlocal electromagnetic currents are given by 

3 



JB-cnAx,y) = Jl /c^a;,4^^(a;i,a;2,a;3)£;g(a:;a:i,a;2,a;3;y), 
^ f 

JB~cm{x,y) = 'Wjdx,J^'^\xi,X2,x:i)E%\x]Xi,X2,x:i]y), 

i—l 
3 

-p-i- /■ dpi f dr ~i J2 Pii^-^i}-ir{x-y) ~ 

EB{x;xi,X2,X3;y) = [[Jj^^Jj^-j^e Eg{pi,p2,P3;r) , 

E^{pi,P2,P3;r) = ^e,,y dT{-w;,i(zi'»i'^"+2(7n*B(--2i)-^«.2K2r"+2<72")$^(-Z2)}, (28) 



Further 



2 = 1 



3 3 

i=l 1=1 

zi = T{wiir + qiY - (1 - T)ql - {wi2r + 92)^ , Z2 = ql - T{wi2r + (72)^ - (1 - t)(?| , 




Wii = 779 : w;'*2 



/ -75(2^2+1^3) 

-j^{2wi+W3) I . (29) 



IV. NUCLEON MASS FUNCTION AND ELECTROMAGNETIC FORM FACTORS 

We start with the calculation of the proton mass function (also called self-energy function) needed for the imple- 
mentation of the compositeness condition. The relevant term in the expansion of the ^-matrix reads 

52 = i^9% JdxJdyp{x){0\T{Jp{x)Jpiy)}\0)p{y) 

= i I dx I dy p{x)Y,p{x - y)p{y) . (30) 



The corresponding two-loop Feynman quark diagram is shown in Fig. 1. In Eq. pop we have introduced the standard 
notation for the proton mass function 

Ep(x -y)^ igl (0|r{Jp(x) Jp(y)}|0) . (31) 
The Fourier-transform of the mass function Tip{x — y) is given by 



Sp(x-y) = J^e-^P'^^-y^^pip), n^ip)^ jd^xe^P^llpix). 



(32) 



We use the same notation Sp for the mass function of the proton in coordinate space and in momentum space. 
Which of the two representations are being used can be read off from the arguments, cif. Sp(a; — y) and Sp(p). The 
matrix element of 52 in Eq. ([30)) between the initial and final proton states (with momenta p and p', respectively) is 
expressed by 

(p'|52b) = i{2T:)U^-^\p-p')up{p)^p{p)uM- (33) 

It is straightforward but nevertheless cumbersome to calculate the proton mass function Sp(a; — y). One uses the 
explicit expression for the interpolating three-quark current given by Eq. ([2]) and the time-ordering of quark fields: 

(0|r{ q%x)qf:{y) }|0) = 5aa'6fr Sf{x - y) ^ 6,,,6ff, j^^^e-^^i-v) Sf{k) , (34) 
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where Sf{x — y) and Sf{k) are the free quark propagators in coordinate and momentum space with 

= 

and where a, a' and /, /' are color and flavor indices, respectively. 
In momentum space the proton mass function is given by 

Sp(p) = I2g% [ / ^!^i,^(_^o)r^^55<i(fci + w.ph^r^tT [Suih - W2p)TASu{k2 ~h+ W3p)rB] , (36) 
J (27r)*iJ [2Tr)H 

where 

^o = 5(fci-fc2)' + ^(fci+fc2)^ (37) 

In order to economize on the notation we introduce a short-hand notation for the two loop momentum integrations 
in ([55)) and in the following formulas. We write 



Note that the integral Eq. (1361) is invariant under a shift of the loop momenta ki — ^ ki + ap where a is an arbitrary 
number and p is the outer momentum. Using this invariance one can obtain various equivalent representations for 
the mass function. In Eq. ([36]) we have chosen a such that the external momentum does not appear in the argument 
of the vertex function. One has W2 = W3 = ^{1 — wi) where wi = md/{md + 2m„). We mention that it is convenient 
to keep m„ 7^ rud in the analytical calculation in order to distinguish the proton from the neutron case. In the end, 
when we do the numerical calculation, we set m„ = rud- 

According to the compositeness condition Eq. ()lip one needs to calculate the derivative of the proton mass function. 
Since the proton is on mass-shell, i.e. u{p)^u{p) = mpfu(j))u{p) and hence p'^u{p)u{p) — mNu(jp)^^u{p), the 
compositeness condition 

S;(p)-^^ = 1 with ^ = m^ (39) 

can be written as 

i9Ep(p) ^ ^.^j^ ^^mN and p^' ^ rriN-f'' ■ (40) 

Here and in the following it is understood that the relations between Green functions are valid when sandwiched 
between spinors. The latter form ([M| is more suitable for our calculation because of its relation to the electromagnetic 
proton vertex function at zero momentum transfer. Using 

d 

— — >5'/(fc + wp) — w Sf{k + wp)^^Sf{k + wp) (41) 



^P^, 



one obtains 



= 123^(($2^(-zo){i«ir^7'^d(fci + wip)rSd{ki + i(;ip)75r^tr [Su{k2 - W2p)T ASu{k2 - fci + W3p)Vb] 

- W2T^J^Sd{ki + Wip)j^T^tT [Suik2 - W2p)Y'Su{k2 - W2p)TASuik2 - fci + W3p)Tb] 

+ W3r^j^Sdiki + wiph^T^'tr [Su{k2 - W2p)TASuik2 ~ k^ + w-ip)i^Su{k2 - fci + W3p)Tb\ })) . (42) 

We now return to the calculation of the electromagnetic vertex of the proton. There are two terms in the relevant ex- 
pansion of the S'-matrix. These are derived i) from the Lagrangian Eq. ([2^ describing the local interaction of the pho- 
ton with the quarks and ii) from the Lagrangian Eq. \21\ describing the nonlocal interaction nucleon+quarks+photon. 
One has 



^3 = i^'gl JdxjdyjdzA^iz)pix){0\T{Jpix) (e„ m(z)7^u(z) rf»7^rf(^)) Jp(2;)}|0)p(y) 

+ i^g% /dx/dy/dz^^(z)p(x)(0|r{Jp(a;) j;_,^(y,z) + j;_,^(x,z)Jp(2/)}|0)p(2/), (43) 
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where the currents Jp{Jp) and Jp-cmiJp-em) ^-re defined by Eqs. ([2]) and (|28|) . respectively. It is important to keep 
track of the signs of the various charges in the calculation. Our choice is to take the electric charges of charged 
particles in units of the proton charge, e.g. Cp = +1, e.^ = +2/3, = —1/3, etc. 

The matrix element Sj, in Eq. (|43p taken between the initial and final proton states with momenta [p) and {p') and 
the photon state with momentum {q — p — p') reads 

{p';q,fJ'\S3\p) = Up{p')T^{p;p',q)up{p) . (44) 

The matrix element T^{p-p'q) is obtained from Eq. ((43l) by the substitutions Af^(z) -> e"^^, p{x) — > ^ and 
p{y) — >■ e~*™. A straightforward calculation gives 

T3^(p;p',9) = z(27r)45W(p-p' - g)A^(p,p') (45) 

where the electromagnetic vertex function A^(p,p') of the proton consists of four pieces represented by the four 
two-loop quark diagrams in Fig. 2: 

• the vertex diagram with the e.m. current attached to the d-quark (Fig. 2a) — 

• the vertex diagram with the e.m. current attached to the u-quark (Fig. 2b) — A^^, 

• two bubble diagrams with the e.m. current attached to the initial proton vertex (Fig. 2c) — ■^p(a)' 
and with the e.m. current attached to the final proton vertex (Fig. 2d) — ^p(b)' 

The four different contributions can be calculated to be 



^pd{P,P') = - 45^((l'Ar(-Zo)*w( - 5(fcl - fc2 + W3qf - i(/ci + ^2 + {2w2 + ^3)?)^) 

X r^7^5d(fci + Wip')rSd{ki + wip' + <7)75r^tr [Su{k2 - W2p')TASu{k2 - k, + W3p')Tb] )) , 

A^'u(p,p') = 16<7ir(($w(-2o)$w( - ^ki -k2- (wi + W2)qf - \{k^ +k2- {wi - W2)qf^ 

X r^7^5rf(fci + i(;ip')7'r^tr [5„(fc2 - i«2p')rA^u(fc2 - fci + W3p')j''Su{k2 - k^ + w^p' + q)TB] )) , 

^p(a)(P'-P') = l2g%(J^^N{-zo)E'^{ki+wip ,-k2 + W2P ,k2-ki+W3p\q) 

X r^7^5d(fci + u;ip')7^r^tr [Su{k2 - W2p)TASu{k2 - fci + W3p')Tb] )) , 

^p{b)(P^P') = 12g^^<^$Ar(-zo)-E^(fci + wip, -A:2 + W2P, fc2 - fci + W3p; -q) 

X rV^<i(fci + '«ip)7'r^tr[5„(fc2-w^2P)rA5„(fc2-fci + ii'3P)rB])). (46) 

The two bubble diagrams Figs. 2c-2d can be seen to be related by 

KibMP')-Kia)(P''P)- (47) 

Gauge invariance requires the validity of the q'^ — Ward identity 

-A^(p,p), (48) 



d^pjp ) 
dptj 



where A^(p,p) is given by the sum of the four contributions Eq. (1461) . The l.h.s. of (1481) has been written down in 
Eq. (l42t . The contribution of the bubble diagrams to the r.h.s. of (|48|) is calculated by using the explicit representation 
of E!^ in Eq. (j28p . For the proton the quark charges are given by ei — and 62 = 63 = e„. One finds 

Kia)(P^P^+Kib^P^P^ - 8{l + 3w,)gl{{k'^^'^{-zonN{-zo) 

X r^^^Sdiki+wip)^^T''tv[Su{k2~W2p)rASu{k2-ki+W3p)rB])). (49) 
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Superficially such terms are not present in Eq. (l42l) since p2|) contains four quark propagators as a result of having 
differentiated the vertex function as compared to the three propagators in P^j) . However, the bubble contributions 
may be rewritten in terms of the vertex diagram contributions at g = 0. This is achieved by using an integration- by- 
parts (IBP) identity where on differentiates the integrand w.r.t. the loop momentum fci. One has 

(( ^{$^(-zo)rV5</(fci + wiph'r^'tT [Suik2 - W2p)rASu{k2 - fci + W3p)Tb] } )) = . (50) 

Upon differentiation and use of the symmetry of the integrand under fc2 — > — ^2 + fci one obtains 

2 (J(Jfc5' ^Ni-zo) $jv(-^o)r^7'^d(fci + wiph^T^'tT [Su{k2 - W2p)TASu{k2 - fci + u'3P)rs] )) 

= (($^(-zo)r^755d(fci + wiph^Sdiki + wip)-f'r''tT [Su{k2 - W2p)rASu{k2 - fci + W3p)Tb] )) 

- ((#^(-zo)r^7'^d(fci + «;ip)7'r^tr [Su{k2 - W2p)rASuik2 - fci + W3Ph''Su{k2 - fci + W3p)rB])) . (51) 

Using the IBP-idcntity and summing up all contributions on the r.h.s. of Eq. (|48l) one finds agreement with the 
l.h.s. of Eq. (|^5)) as given by Eq. (H^ . One has thus proven the validity of the Ward identity (|48p (recall that 

W2 = W3 = 1/2(1 - Wi)). 

A further technical remark is in order concerning the above check of the Ward identity Eq. (|48p . The proof made 
use of an IBP identity assuming the vanishing of the pertinent surface term. However, in the confinement ansatz with 
the accompanying IR cutoff the requisite surface terms no longer vanish. As it turns out one can avoid the use of IBP 
identities in the proof of the Ward identity by astutely shifting the loop momentum fci in the mass function. The 
appropriate shift is fci fci + (w2 + — |)p and fc2 — fc2 + W2P- One then obtains 

Sp(p) = 12g^ J 1 -01-3%{-z,)T^j'S4ki - ip)7'r^tr [5„(fc2)rA5«(fc2 - fci + |p)rs] , (52) 

where 

Zi = i(fci - fc2 + iw3 - l)pf + i(fci + fc-2 + (2W2 +W3~ l)pf . (53) 

After the shift of the loop momentum the argument of the vertex function $Ar(— zi) now depends on the external 
momentum p. When differentiating w.r.t. the external momentum p a new term will appear caused by the derivative 
of the vertex function $jv in addition to the terms originating from the derivatives of the quark propagators. After 
shifting back the loop momenta zi — zq and some algebraic juggling one finds that the derivative of the mass function 
coincides analytically with the expression for the electromagnetic vertex function A^{p,p) given by the sum of the 
contributions of the triangle diagrams and the bubble diagrams Eq. (HHl) . We emphasize that in this derivation 
we did not have to make use of an IBP identity to prove the Ward identity. 
Hereafter we will use the compositeness condition = in the form 

A^d(p,p)+A^„(p,p)-fA^(^^(p,p)+A^(,^(p,p) =7^' with {fi = mN) and = tojv7^ (54) 

in order to determine the coupling constant g^- This allows one to provide the correct normalization of the charged 
proton form factor within the confinement scenario. 

Another useful check is to reproduce the generalized Ward-Takahashi identity 

q^A^'(p,p') = Ep(p)-Ep(p'). (55) 

We shall not elaborate on this proof which is straightforward by using suitable shifts of the loop variables. 

Let us briefiy describe another check on the gauge invariance of our calculation. Without gauge invariance there 
are three independent Lorentz structures in the electromagnetic proton vertex which can be chosen to be 

A^^ip,p') ^ r Ffiq^) - ^ F!{q') + Ff,^{q^) , (56) 

where cr'"? = '^{l^l'^ ~ l'^l^)<lv- The form factor Ff^Q{q^) characterizes the non-gauge invariant piece and must 
therefore vanish for any q^ in a calculation which respects gauge invariance. For the four contributions of Fig. 2a-2d 
we found that 

i^^c,(<Z^)^0, ^^^^Jg2)^0, Fl^^,^Aq-)^^Fl^^,y) ^q\ (57) 
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The gauge variant contributions of the two vertex diagrams are zero while they vanish for the sum of the two bubble 
diagrams. 

Before discussing the e.m. properties of the neutron we would like to comment on a potential conflict between gauge 
invariance and our confinement ansatz. In general the IR cutoff used in Eq. (jl9p can destroy the gauge invariance 
as any cutoff can do. One can, however, show that in some special cases gauge invariance remains unimpaired when 
implementing confinement through an IR cutoff. For example, in Appendix B of [HI we have shown that the p — 7 
transition amplitude is gauge invariant off mass-shell even in the presence of an IR cutoff. The crucial point of the 
proof was that we were able to show that the integrand of the t-integration itself was gauge invariant. In the case of 
the electromagnetic form factor of the proton one finds again that the integrand of the t-integration is gauge invariant 
by itself due to a symmetry property of the integrand in the space of the Schwinger a-parameters. However, if the 
proof of gauge invariance requires an integration by parts in the space of momenta which becomes translated into 
an integration by parts over the t-parameter gauge invariance will be spoiled by the surface term due to the upper 
integration limit 1/A^. In order to keep gauge invariance one can proceed as follows. First, by using the properties of 
the relevant integrals over the loop momenta one needs to specify a gauge invariant part of the full amplitude. Then 
one employs our confinement ansatz for the gauge invariant parts of the amplitudes. Such an approach was used to 
verify the validity of the Ward identity when connecting the derivative of the mass function and the electromagnetic 
vertex function in the presence of an IR-cutoff. 

The electromagnetic vertex function of the neutron is obtained from that of the proton by replacing m„ -f-^- rrid, -f-^- 
Cd and rup — >■ m„. F-j^ {q^) and (q^) arc the Dirac and Pauli nuclcon form factors which are normalized to the electric 
charge cn and anomalous magnetic moment {kN is given in units of the nuclear magneton e/2mp), respectively, 
i.e. one has F^' (0) = ejv and i^2^(0) = kN- In particular, one can analytically check by using the IBP identity that 
the Dirac form factor of the neutron is equal to zero at = 0. 

The nucleon magnetic moments fiN = -P"i^(0) + -^2^(0) are known experimentally with high accuracy [20| 

Ai^'^P* = 2.79 Ai^'^P' = -1.91 . (58) 

We will use these values to fit the value of the nucleon size parameter. The other model parameters are taken from 
the fit to mesonic transitions done in lal: 



rrib 



(59) 



0.235 0.424 2.16 5.09 0.181 GeV 

We obtain 



vector current ^ Aat = 0.36 GeV Hp ^ 2.79 /x„ = -1.70, (60) 

tensor current Ajv = 0.61 GeV fip — 2.79 /x„ = — 1.69. (61) 

It is convenient to introduce the Sachs electromagnetic form factors of nucleons 



G^iq') = F,^iq^) + /^Fi'iq^), 



GZil') = F-{q^) + F-{q^). (62) 
The slopes of these form factors are related to the well-known electromagnetic radii of nucleons: 



2\N _ (.""^N 



dG%{q^)_ 

N f„2\ 



dq^ 

6 dGUq^ 



^ ''^ Gm dq 



(63) 
(64) 



We would like to emphasize that reproducing data on the neutron charge radius (?'|;)" is a nontrivial task (see e.g. 
discussion in Ref . [lOj ) . As well-known the naive nonrelativistic quark model based on SU(6) spin-flavor symmetry 
implies = 0. The dynamical breaking of the SU(6) symmetry based on the inclusion of the quark spin-spin 

interaction generates a nonvanishing value of From this point of view the dominant contribution to the {r%)" 

comes from the Pauli term: 



12 



The experimental data on the nucleon Sachs form factors in the space- Hke region = — > can be approximately 
described by the dipole approximation 

1 + ^ M« ^ M« ^ (l-gVO.TlGeV^)' 

According to present data the dipole approximation works well up to 1 GeV^ (with an accm-acy of up to 25%). For 
higher values of the deviation of the nucleon form factors from the dipole approximation becomes more pronounced. 
In particular, the best description of magnetic moments, electromagnetic radii and form factors is achieved when we 
consider a superposition of the V- and T-currents of nucleons according to Eq. ([3]) with x = 0.8. For the size 
parameter of the nucleon we take Ajv =0.5 GeV. 

In Table I we present the results for the magnetic moments and electromagnetic radii for this set of model parameters. 
In Fig. 3 we present our results for the dependence of electromagnetic form factors in the region G [0, 1] GeV^. 
Fig. 3 also shows plots of the dipole approximation to the form factors. The agreement of our results with the dipole 
approximation is satisfactory. Inclusion of chiral corrections as, for example, developed and discussed in 21] may lead 
to a further improvement in the low description. 



V. A-TYPE MASS FUNCTION AND ELECTROMAGNETIC VERTEX 

In a future publication we plan to study the rare baryon decays Af, — > A^ £^£^ in the context of the covariant quark 
model |9] . It is the purpose of this section to provide the necessary material that allows for a covariant quark model 
description of the A — {Q[ud])-type baryons composed of a (s,c, &) quark Q and a light diquark-like state [ud] with 
spin and isospin zero. 

In general, for the A-type baryons one can construct three types of currents without derivatives — pseudoscalar 
J^, scalar and axial- vector (see, Refs. P,^,!!!!!!): 

JrQ[„,] = e'^^'^^'^^Q'^'u^'Cj^d''', 

•^Aqm = e"^"^"^7^Q°^""^C757Mt^°'- (67) 

There are only two independent linear combinations of the above three currents given by — (2J^ — 2J^ + J^)/3 
and J"^ = + . The symbol [ud] denotes antisymmetrization of both flavor and spin indices w.r.t. the light 
quarks u and d. We will consider three flavor types of the A-baryons: Ag[u(i], A+[ud] and AjJ[ud]. In Ref. [3|] we 
have shown that, in the nonrelativistic limit, the and J'^ interpolating currents of the Aq[^,j;] baryons become 
degenerate and attain the (same) correct nonrelativistic limit (in the case of single-heavy baryons this limit coincides 
with the heavy quark limit), while the current vanishes in the nonrelativistic limit. On the other hand, the 
and interpolating currents of A-type baryons become degenerate with their SU(Nf )-symmetric currents in the 
nonrelativistic limit. In Ref. [1^ we have shown that in case of the heavy-to-light baryon transition A+ A^e'^Ve 
the use of a SU(3) symmetric current for the A° hyperon is essential in order to describe data on r(A+ A'^e+i^e) 
(see also discussion in Refs. [l^, [l^). Therefore, in the following we restrict ourselves to the simplest pseudoscalar 
current. The nonlocal interpolating three-quark current is written as 

Ja{x) = dxi dx2 dx3FA{x;xi,X2,X3) J^^\xi,X2,X3) , (68) 



'3q 

4^\xuX2,X3) = ^e'^'"'''' QaA^i) {u''Hx2)C-i'>d-'{x^)-d''-{x^)C-i^u'^Hx2)) 
= e"^"^"-^ Q^'ixi) ii"^(a;2) df^'ixz) , 

Jk{x) = JdxiJdx2jdx3F/^{x;xi,X2,X3) jI^^\xi,X2,X3) , 
4^\xi,X2,X3) = e^'^-^-d'''{x3)^''Cu'''{x2)-Q^'{xi), 

where Q = s,c,b. 

The calculation of the A-type mass function and the electromagnetic vertex proceeds in the same way as in the 
nucleon case. The matrix elements in momentum space read 

Sa(p) = 6gl(^(^^l{-zo)SQ{h + wip)tT [5„(fc2 - W2Ph^Sdik2 -h+ wsp^'] )) , (69) 
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where we use the same short-hand notation << ... >> for the two- fold loop- momentum integration as before (see 
Eqs. dM])). The variable zq is defined in ([57|. 

The various contributions to the electromagnetic vertex are given by 

KIq{p,p') = 6eQgi<J^I>A(-zo)*A( - 5(^1 - ^2 + w^qf -\{ki + k2 + (2w2 + u-s)?)^) 

X SQ{ki + wip)^^SQ{ki + wip + q)tr [S'„(fc2 - W2p')7^ 5*^(^2 - ki + W3p')j^] , 

KuiP^P') = - 6e„5i^<J$A(-.2o)^A( - - A:2 + w^qf - ^(/ci + fc2 - (2i(7i + u-g)?)") 

X S'Q(fci + u;ip')tr [S'„(fc2 - W2J3' - q)l^Su{k2 ~ W2p')^^'' Sd{k2 - fci + wspOt^] )) i 

KdiP'P') = 6ed5A((^A(-2o)$A( - 5(^:1 - fca - K W2)qf - ^{ki + ka ~ {wi - ^2)9)^) 

X Sqiki + wip')tv [Su{k2 - W2p'h^Sd{k2 - fci + W3p')-f'' Sd(k2 - fci + w^p' + 9)7^] J) 

^A(a)(P'P') = 6gi<^(^$A(-zo)£'A(fci +u;ip',-fc2 +-u;2p',fc2 - fci + wsp';?) 
X Sgiki + wip')tr [Su{k2 - W2p'H^Sd{k2 ~ ki + wsp')^ 



^A{b)(P'P') = GgA{{^A{-ZQ)E'^{ki+wip,-k2 + W2P,k2-ki+W3p;~q) 

X Sqiki + wip)ii [Su{k2 - W2p)i^Sd{k2 - fci + W3P)7^] )) • (70) 

One now has three e.m. vertex contributions because there are three different quarks in the Ag state. The function 
E'^{'ri,r2.r'i\r) has been defined in Eq. ([28)) . The variables qi = X]i=i ^ii''^ ^.nd 52 = Si=i ^i2?'i m £;j^(ri, r2.r3; r) 
can be seen to be related to the loop momenta by 

<Zi - 75(fci -fc2), 92 = -75(fci + fc2) (71) 

for both bubble diagrams. By using Eq. ([7T]) one finds the q — Q relations 

A5((^)(p,p)-|-A^(,)(p,p) = -8gi(((Qifcf -f Q2fc2^)$k(-^o)§A(-^o) 

X Sqiki + wip)iv [5„(fc2 - W2p)^^Sd{k2 - fci + ^3^)7^] 



Qi = ei(w2 + 2^3) — 62(101 — W3) — 63(2^1 -I- u;2) , 

Q2 ~ ei(w2 — W3) — 62(^1 2W3) -f 63(^1 + 2^2) • (72) 

where the subscripts on the charges refer to the flavors of the three quarks: "i = 1" — !• "s,c, 6", "i = 2" — ;> "u" and 
" i = 3" — >■ " d" . Next we will use an IBP-identity to write 

d 
-Sfcf 

One finds 



((^{*i(-2o)^Q(fci+i^ip)tr[5„(fc2-i«2P)7'5rf(fc2-fci+u;3P)7'] })) = 0, (« = 1,2). (73) 



where 



fc^'^o)) = i(((Ai' + 2A^ + A^'))), (74) 

^0 = *a(--zo) *A(-^o)5'Q(fci wip)tr [S'„(fc2 - W2p)7^5'd(fc2 - fci + 103^)7^] , 

^1 = $A(-2o)S'Q(fci -I- wip)7^5Q(fci -I- u>ip)tr [5'„(fc2 - u;2p)7^'5'd(fc2 - fci + W3P)7^] , 

^2 = $i(-^:o)5'Q(fci+wip)tr[S'„(fc2-W2p)7''5'«(fc2-W2p)7^S'<i(fc2-fci-Hw3P)7^] , 

^3 = ^\{-zo)SQ{ki + wip)tr [^^(fca - W2p)7^'S'd(fc2 - fci -I- z«3P)7'''S'<i(fc2 - fci + W3P)7^] • (75) 
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Using these identities and collecting all pieces together, one has 

Aa(p,p) = (eg + e„ + e<j) ^^^^ , ^ = toa . (76) 

As was discussed above, this Ward identity allows one to use the compositeness condition = written in the form 

KiP^p)=r, ?^ = ™a, (77) 

where we take eg = for the present discussion. Again we have checked analytically that, on the A-type baryon 
mass shell, the vertex diagrams are gauge invariant by themselves and the non-gauge invariant parts coming from the 
bubble diagrams corresponding to Fig. 2(c) and 2(d) cancel each other before i-integration. The standard definition 
of the electromagnetic form factors is 

A>;^{p,p')=j^F,{q^)^—.F2{q'), (78) 
where a^'^ — 1(7^7"^ — l'^l^)qv The magnetic moment of the A-type baryon is defined by 

MA = (^^i(0)-f F2(0)) (79) 

ZTOA 

In terms of the nuclear magneton (n.m.) 2^ the A-type baryon magnetic moment the A-hyperon magnetic moment 
is given by 

MA - (Fi(0) + i^2(0)) ^, (80) 
mA 

where mp is the proton mass. 

In the present paper we shall only make a rather cursory investigation into the possible values of the size parameters 
of the A = (Q[uc?])-type baryons. A more detailed investigation will be left to our future publication 9] where we 
will include information on the charged current transitions A^ — > A^ and Ac A^ to specify the values of the size 
parameters of the A = {Q[ud])-iyY>^ baryons. 

Let us assume for the moment that the size parameters are the same for all A-type baryons. One then has 



MA, = 


-0.73, 


MA, = 


+0.36, 


MA, = 


-0.06, 


MA, = 


-0.68, 


MA, = - 


-fO.40, 


MA, = 


-0.06, 


MA, = 


-0.61 , 


MA, = 


+0.44, 


MA, = 


-0.07. 



(81) 

The magnetic moment of the A^ has to be compared with the experimental value listed in (20j 

MA, = -0.613 ± 0.004. (82) 



Eq. (|81[) shows that the value of the magnetic moment of the As is quite stable against variations of its size parameter. 
There is no experimental information on the magnetic moments of the Ah and Ac. 

The calculation of the form factors in our approach is automated by the use of FORM and FORTRAN packages 
written for this purpose. In order to be able to compare with our earlier unconfined calculations we have written two 
versions for the confined and the unconfined versions of the covariant quark model. 



VI. SUMMARY AND CONCLUSIONS 



We have extended our previous formulation of the confined covariant quark model for mesons and tetraquark states 
to the baryon sector. We have discussed in some detail various calculational aspects of the two-loop baryon problem 
such as the evaluation of the baryon mass operator and its derivative, the implementation of confinement in the 
two-loop context, the calculation of electromagnetic current-induced transition matrix elements and the analytical 
verification of the pertinent Ward and Ward-Takahashi identities associated with the electromagnetic matrix elements. 

In our numerical work we have used the same values of the constituent quark masses and infrared cutoff as had 
been obtained before in the meson sector by a fit to various mesonic transition matrix elements. In this way the 
number of model parameters were kept to a minimum. 
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Using two parameters we have calculated the nucleon magnetic moments and charge radii as well as the electro- 
magnetic form factors at low momentum transfers. An extension of our work to the N — A(1236) transition can be 
done along the lines described in [26l |. 

We have also discussed light and heavy A — (Q[uc?])-type baryons. In particular we obtained a value for the size 
parameter of the A^ by a fit to its experimentally known magnetic moment. By determining the properties of the 
A = {Q[ud])-type baryons we have laid the groundwork for a calculation of the rare decays of the Ab-baryon (such as 
Ah — As£~^i~) within the framework of the covariant quark model. 
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TABLE I: Electromagnetic properties of nucleons 



Quantity 


Our results 


Data [20] 


fip (in n.m.) 


2.96 


2.793 


^j,„ (in n.m.) 


-1.83 


-1.913 


(fm) 


0.805 


0.8768 ± 0.0069 


(r|)" (fm^) 


-0.121 


-0.1161 ± 0.0022 


rlj (fm) 


0.688 


0.777 ± 0.013 ± 0.010 


rli (fm) 


0.685 


n sRo+O-009 
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FIG. 2; Electromagnetic vertex function of the proton: (a) vertex diagram with the e.m. current attached to d-quark; (b) 
vertex diagram with the e.m. current attached to u-quark; (c) bubble diagram with the e.m. current attached to the initial 
state vertex; (d) the bubble diagram with e.m. current attached to the final state vertex. 
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FIG. 3: Sachs nucleon form factors in comparions with the dipole representation in the space-hke region Q < 1 GeV^ 



